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Abstract 

We consider the neutral triple gauge boson production process in the context 
of large extra dimensions (LED) model including the Kaluza-Klein (KK) excited 
gravitons at future linear colliders, say ILC(CLIC). We consider 777, 77Z, -yZZ 
_ and ZZZ production processes, and analyse their impacts on both the total cross 

i-^ 1 section and some key distributions. These processes are important for new physics 

searches at linear colliders. Our results show that KK graviton exchange has the 
most significant effect on e~e + —> jZZ among the four processes with relatively 
small Ms, while it has the largest effect on e~e + — > 777 with larger Ms- By using 
the neutral triple gauge boson production we could set the discovery limit on the 
fundamental Plank scale M5 up to around 6-9 TeV for 5 = 4 at the 3 TeV CLIC. 
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1 Introduction 



The hierarchy problem of the standard model (SM) strongly suggests new physics at 
TeV scale, and the idea that there exists extra dimensions (ED) which first proposed by 
Arkani-Hamed, Dimopoulos, and Dvali[T] might provide a solution to this problem. They 
proposed a scenario in which the SM field is constrained to the common 3+1 space-time 
dimensions ("brane"), while gravity is free to propagate throughout a larger multidimen- 
sional space D = 5 + 4 ("bulk"). The picture of a massless graviton propagating in D 
dimensions is equal to the picture that numerous massive Kaluza-Klein (KK) gravitons 
propagate in 4 dimensions. The fundamental Planck scale Ms is related to the Plank 
mass scale M Pl = G^ 1/2 = 1.22 x 10 19 GeV according to the formula M 2 Pl = 8irM s s +2 R 5 , 
where R and 5 are the size and number of the extra dimensions, respectively. If R is large 
enough to make Ms on the order of the electroweak symmetry breaking scale (~ 1 TeV), 
the hierarchy problem will be naturally solved, so this extra dimension model is called the 
large extra dimension model (LED) or the ADD model. Postulating Ms to be 1 TeV, we 
get R ~ 10 13 cm for 5 — 1, which is obviously ruled out since it would modify Newton's 
law of gravity at solar-system distances; and we get R ~ 1 mm for 6 = 2, which is also 
ruled out by the torsion-balance experiments [2]. When 5 > 3, where R < 1 nm, it is 
possible to detect graviton signal at high energy colliders. 

At colliders, exchange of virtual KK graviton or emission of a real KK mode could 
give rise to interesting phenomeno logical signals at TeV scale [31 H]. Virtual effects of 
KK modes could lead to the enhancement of the cross section of pair productions in 
processes, for example, di-lepton, di-gauge boson (77, ZZ, W + W~), dijet, ti pair, HH 
pair [5], El [3, El [9], QUI E] etc. The real emission of a KK mode could lead to large missing 
Et signals viz. mono jet, mono gauge boson[3l HI [T2], [13] etc. The CMS Collaboration 
has performed a lot of search for LED on different final states at yfs = 7 TeV[T4"l IT5| [16], 
and they set the most stringent lower limits to date to be 2.5 TeV < Ms < 3.8 TeV by 
combining the diphoton, dimuon and dielectron channels. 

Studies for LED have been extended to three body final state processes in recent years. 
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Triple gauge bosons productions in the SM are important because they involve 3-point 
and 4-point gauge couplings in the contributing diagrams, which allow for restrictive tests 
of triple and quartic vector boson coupling. And also they might contribute backgrounds 
to new physics beyond the SM. Furthermore they are sensitive to new physics. This 
kind of processes have been studied at LO[TH HSJ US] and NLOJ2D1 [2D 1221 1231 [21 122] in 
the SM, and virtual graviton exchange effects to these processes within the LED model 
at LHC are studied recently [25]- Linear colliders have more advantage in testing extra 
dimensions than LHC for the following reasons. First, even though the LHC has much 
higher center-of-mass (c.m.s.) energy than linear colliders, the theoretical amplitude at 
LHC is hampered by the unitary constraint a/I < Ms. Second, linear colliders have cleaner 
environment than LHC, so it's much easier to select the ED signals. The capabilities of 
the planned International Linear Collider (ILC) and Compact Linear Collider (CLIC) for 
precision Higgs studies are well documented \27\ |2"8] . They will also provide opportunities 
for the search for new physics beyond SM. So in this paper we consider the triple gauge 
bosons production at ILC and CLIC e~e + — > VVV within the LED model, where we 
restrict V to be neutral gauge boson (V = "f,Z). The following four final states are the 
subject of this analysis: (i) 777 (ii) 77Z (iii) jZZ (iv) ZZZ. The case where V = 
is in prepare and will be part of a different paper. 

This paper is organized as follows: in section 2 we present the analytical calculation 
of the processes mentioned above with a brief introduction to the LED model, section 3 
is arranged to present the numerical results of our studies, and finally we summarize the 
results in the last section. 

2 Theoretical Framework 

In this section we give the analytical calculations of the process e~e + — > VVV with 
V = 7, Z at linear colliders in the LED model. In our calculation we use the de Donder 
gauge. The relevant Feynman rules involving graviton in the LED model can be found in 
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Ref.jl]. We denote the process as: 

e'ivx) + e + (p 2 ) V{p 3 ) + V(p,) + V(p 5 ) (1) 



where pi, p 2 and p$, p±, p§ represent the momenta of the incoming and outgoing particles 
respectively. 




(a) (b) (c) (d) (e) 

Figure 1: Feynman diagrams for e~e + — > VVV process with V = 7, Z. (a,b) are the 
diagrams in the SM, and (c,d,e) are the diagrams in the LED model. 



In Figfj] we display the Feynman diagrams for this process in both the SM and LED 
model, among which (a) and (b) are SM diagrams and (c)~(e) are LED diagrams. We have 
neglected the Higgs coupling to electrons because the Yukawa coupling is proportional 
to the fermion mass. But the Higgs coupling to Z bosons, which appears in e~e + — > 
ZZZ process (Figfjjb)), can't be neglected because of its large contribution, e.g., with 
y/s to be 300 - 3000 GeV, the cross section for FigJT](b) is about 26% - 9% of the total 
SM cross section. After considering all possible permutations, we have 12 SM diagrams 
for the ZZZ production process and 9 SM diagrams for the other three processes, and we 
have 12 LED diagrams for the 777 and ZZZ processes and 4 LED diagrams for 77Z and 
jZZ processes. 

In our calculation we consider both the spin-0 and spin-2 KK mode exchange effect. 
The spin-0 states only couple through the dilaton mode, which have none contribution 
to 777 and 77Z processes and could contribute to ^ZZ and ZZZ production processes 
through couplings to massive gauge bosons. However the cross sections coming from the 
dilaton mode are so small that can be neglected, e.g., they are at most about 10 -12 and 
10 -8 times of the total cross sections for 7ZZ and ZZZ production processes, respectively. 
So we focus our study on the spin-2 component of the KK states. 
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The couplings between gravitons and SM particles are proportional to a constant 



named gravitational coupling k = \/16tiGn, which can be expressed in terms of the 
fundamental Plank scale Ms and the size of the compactified space R by 

k 2 R & = 87r(47r) 5 / 2 r(5/2)M~ {5+2) (2) 

In practical experiments, the contributions of the different Kaluza-Klein modes have 
to be summed up, so the propagator is proportional to if (sy — where = (pi+Pj) 2 
and rrifi is the mass of the KK state n. Thus, when the effects of all the KK states 
are taken together, the amplitude is proportional to — _ % = D(s). If 5 > 2 this 
summation is formally divergent as becomes large. We assume that the distribution 
has a ultraviolet cutoff at ~ Ms, where the underlying theory becomes manifest. 
Then D(s) can be expressed as: 



D ^ = i^w s )S ~ 2l * +2iI{Ms/v ~ s)] - (3) 



The imaginary part l(A/y/s) is from the summation over the many non-resonant KK 
states and its expression can be found in Ref. [1]. Finally the KK graviton propagator 
after summing over the KK states is: 

Using the Feynman rules in the LED model and the propagator given by Eq.(jH), we 
can get the amplitudes for the virtual KK graviton exchange diagrams in FigJTJ The 
total amplitude can be obtained by adding these LED amplitudes together with the SM 
ones. The total cross section can be expressed as the integration over the phase space of 
three-body final state: 



(2tt^ 4 

0~tot 



/ d<S> 3 J2\M tot \ 2 . (5) 



4|pi 

where ^2 represents the summation over the spins of final particles and the average over 
the spins of initial particles. The phase-space element c?$ 3 is defined by 



3 Numerical Results 

3.1 Input parameters and kinematical cuts 

We use FeynArts and FormCalc package [23 ED] to generate and reduce the amplitudes 
and then implement numerical calculation. We use BASES [31] to perform the phase space 
integration and CERN library to display the distributions. The SM parameters are taken 
as follows [32]: 

m z = 91.1876 GeV, m w = 80.399 GeV, m e = 0.511 MeV, 
a{m%) = 1/127.934, m H = 125 GeV[33j E3|. 

We take the cuts on final particles as: 

vf > 15 GeV, r/ 7 ' z < 2.5, # 77 ^ 0.4 (7) 

where R is defined as R = a/ (A0) 2 + (A?]) 2 , with A0 and Arj denoting the separation 
between the two particles in azimuthal angle and pseudo-rapidity respectively. 

3.2 Total Cross sections 

In Figj2]we present the cross sections for e~e + — > 777, 77Z, 7ZZ and ZZZ processes as 
the functions of yfs for 5 = 3 with different values of Ms- The solid lines are the SM 
results. The dashed, dotted and dot-dashed lines are corresponding to the cross sections 
in the LED model for Ms = 3.5 TeV, 4.5 TeV and 5.5 TeV respectively. When yfs is 
less than about 1 TeV, the curves for the total cross sections including the LED effect 
seem to be overlapped with that in the SM, then the LED effect becomes significant 
with the increment of \fs. When Ms = 3.5 TeV, the e~e + — > ^yZZ process has the 
most significant LED effect among the four processes considered in this paper, while the 
e~e + — > 77 Z process has the least. When y^s = 500 GeV, the cross sections for 777 and 
77Z production processes are comparable, which are 91 fb and 78 fb respectively, because 
they have comparable phase-space. While the cross sections for 7ZZ and ZZZ processes 
are much smaller due to the less phase-space, which are 18 fb and 1 fb, respectively. With 
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yfs increase to be 3 TeV, the cross section for ^ZZ are enhanced to 347 fb, which is even 
larger than the 777 process (323 fb), and the cross sections for ZZZ process is enhanced 
to 34 fb, which is comparable to the 77Z process (36 fb). With larger Ms value (5.5 
TeV), the LED contribution to 777 production process will exceed ^ZZ , that's why 
e~e + —> 777 process puts the highest limits on M s , as we will see later. 
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Figure 2: The cross sections for the process e e + — > VVV with V = 7, Z in the SM and 
LED model as the function of with M s = 3.5, 4.5, 5.5 TeV and 5 = 3. 



In Figj3] we present the dependence of the cross section on energy scale Ms with 
y/s = 1 TeV, 2 TeV and 3 TeV respectively. In each figure of FigSJ, we present the curves 
for the cross sections with the extra dimension 5 value being 3, 4, 5 and 6 separately. The 
solid straight lines, which are independent of Ms, are the SM results, and the dashed, 
dotted, dash-dotted and dash-dot-dotted lines are the cross sections for 5=3, 4, 5 and 6 
respectively. It's clear that for a given value of 5, the cross section decreases rapidly with 



the increment of Mg, and finally approaches to its corresponding SM result. We can see 
again that the virtual KK graviton exchange contribution decreases with the increment 
of the 5 value. The LED effect on the cross sections with y/s = 1 TeV is too small to be 
detected, especially for e~e + — > 77Z process, which is coincidence with FigfSJ If we got 
high enough c.m.s energy, say 3 TeV, the cross sections would be very significant when 
Ms is not very large. Even with Ms = 6 TeV, the cross sections are still several times of 
the SM ones. 
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Figure 3: The cross sections for the process e e + — > VVV in the SM and LED model as 
the function of Ms with \/s = 1, 2, 3 TeV and 5 = 3, 4, 5, 6. 



3.3 Distributions 

The distributions of the Gauge boson pair invariance mass Myy (VV = 77, ZZ) and the 
Gauge boson transverse momentum as well as their rapidity y v at the 3 TeV CILC, are 



S 



shown in Fig.4-6. The results are for Ms = 6.5 TeV at the fixed value 4 for the number 
of extra dimensions and obtained by taking the input parameters mentioned above. 



e e + — > 777 

Before selecting our event samples for triple 7 production, we order the photons on the 
basis of their transverse momentum i.e., p^ > > p^. For e~e + — > 777, we are 
interested in the p\ and y 1 distribution which are displayed in the left and right panel in 
Fig.4, respectively. The solid, dashed and dotted lines refer to 71, 72 and 73, respectively. 
In high p^ and p^ region, the LED effect dominant the total (SM+LED) distribution, 
because more KK modes contribute with the increase of pt- Difference can be found for 
the production, although it's still enhanced by the LED effects, it's low px region is 
dominant while in high px region it becomes much smaller. Rapidity distribution of the 
related photon has been shown in the right panel in Fig.4. As we can see, the rapidity 
distributions in the LED model show significantly peaks around y = 0, which implies the 
large contributions at high p T 1,2 region. Compare with the 71 and 72 distribution, the y 
distribution of 73 seems much flatter. 
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Figure 4: The transverse momentum(p^) and Rapidity(y) distribution of photons for the 
process e~e + — > 777, on the basis of their transverse momentum p T x > p T 2 > pl^ with 
M s = 6.5TeV, y/s =3 TeV and 5 = 4. The solid, dashed and dotted lines refer to p^iy 11 ), 
Pt 2 ^ 2 ) andp]?(y 73 ). 
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e-e + -> ZZZ 



Similar to the e~e + — > 777 production, triple Z bosons final particles are classified in 
such a way that p^ 1 > p? 2 > p^ 3 - Similar conclusion can be found for the e~e + — > ZZZ 
production. It's not strange that the signal of triple 7 signal is larger than the triple Z 
production since the three Z bosons suppress the phase space integration extremely, so 
that the total cross sections as well as the distributions become smaller as can be seen 
in Fig. 4 and Fig. 5, the peak is around 0.0006 fb/GeV for p^ 1 compared to 0.01 fb/GeV 
for p^ in the high px region. For the y distributions, the peaks for the ZZZ production 
are narrower than the 777 distributions, however, the conclusion is the same that the 
rapidity distributions in the LED model show significant peaks around y = 0. 
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Figure 5: The transverse moment um(j>r) and Rapidity(y) distribution of Z bosons for the 
process e~e + — > ZZZ, on the basis of their transverse momentum p^ 1 > p? 2 > Vr' with 
Ms = 6.5 TeV, \fs = 3 TeV and 5 = 4. The solid, dashed and dotted lines refer to 
PT(y Zl ),PT 2 (y Z2 ) and p^(y z *). 



e e + — > "ffZ and e e + — > "fZZ 

Now let's see the distributions for the e~e + —> 77Z and e~e + — > ^yZZ productions. The 
photon pair decay of the KK graviton is one of the clean decay modes, so the distribution 
of the invariant mass of the photon pair (M 77 ) is a useful observable for e~e + — > 77Z . 
An obvious enhancement on the tail of this distribution makes such region of extreme 
interest. Typically, we find that the KK modes dominate over SM contribution for larger 
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values of invariant masses (say above 1 TeV for a given set of Ms and S values, here we 
give Ms = 4.5 TeV and 5 = 4) of photon pairs indicating the observable nature of the 
signal, see the two dotted line in Fig.6. The upper and lower ones refer to the SM predict 
and SM+LED effects. For the process e~e + — > jZZ, it is similar to 77Z production 
process, and in this case, the invariant mass of Z boson pair(Mzz) is a useful observable. 
We thus display it in Fig.6, see the solid lines. These two solid line branch at about 
the invariance mass 1 TeV, and the upper and lower ones present the SM and SM+LED 
effects, respectively. 
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Figure 6: Invariant mass distribution of M^zz) for e e + — \ r y r yZ(e e + —> ^ZZ) for 
M s = 4.5 TeV, y/s = 3 TeV and 5 = 4. 

It is clear that if the deviation of the cross section from the SM prediction is large 
enough, the LED effects can be found. We assume that the LED effects can and cannot 
be observed, only if 



Ao- = \a tot - a SM \ > = 5cr 



and 



a i i ^ 3y/£a tot _ 
l\a = \a tot - a S M\ < = 6a 



(9) 



Our final results show that by using the 777 production we can set the discovery limit 
on the fundamental Plank scale Ms up to 3.1-9.9 TeV, depending on the extra dimension 
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1 TeV 


2 TeV 


3 TeV 






5cr 3a 


5a 3a 


5a 3a 


e~e 4 




3599 4093 


6371 7187 


8906 10023 


e~e 4 


zzz 


2530 2858 


4711 5312 


6730 7584 


e~e 4 


' —> 77Z 


2052 2289 


4007 4458 


5823 6500 


e~e 4 


- -»■ 7ZZ 


3406 3835 


6027 6757 


8412 9421 



Table 1: The discovery (Aa > 5a) and exclusion (Aa < 3a) LED model 
fundamental scale (Ms) values for the e~e + — > VVV processes at the y/s = 1, 
2, 3 TeV ILC(CLIC). C = 300 fb~\ 5 = 4. 

5 C [3,6], with the luminosity 300 fb _1 and the colliding energy 1-3 TeV. For the other 
three final states 77Z , ^ZZ and ZZZ , the limits are 1.8-6.3 TeV, 2.9-9.3 TeV and 
2.2-7.4 TeV, respectively. To do a more detailed description, in Table 1, we present the 
5a discovery and 3a exclusion fundamental scale Ms values at the ILC/CLIC with the 
luminosity 300 fb -1 for 5 = 4. It shows that compared to the other three channels, 777 
can set the discovery limit bounds much higher, up to 8.9 TeV. The phenomenology of the 
neutral triple gauge boson production at the near future is much richer at linear colliders, 
though its production cannot give compete limits as, for example, dilepton production 
gives, it's still very interesting and important. 





C (fb- 1 ) 
5a 


3a 


PP ~^ 111 


960 


1500 


pp — > zzz 


21 


21 


PP ~^ 11% 


20 


17 


PP ~^ iZZ 


170 


140 



Table 2: Integrated luminosity needed at 14 TeV LHC to accomplish the dis- 
covery and exclusion bounds at a 1 TeV linear collider, which are listed in the 
first two columns in Table 1, using the neutral triple gauge boson production 
processes, with 5 = 4. 



To make a comparison with Ref . |26] . we repeat the pp — > VVV (V = j,Z) process 
at LHC, using the same parameters and cuts with Ref. [26], and find that our results are 
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in good agreement with theirs. In Table 2 we list the integrated luminosity the 14 TeV 
LHC needed to accomplish the discovery and exclusion bounds at a 1 TeV LC (the first 
two column data listed in Table 1), by using the corresponding VVV production channels 
with extra dimensions 5 = 4. The table shows that with years of collection of data, LHC 
could accomplish the discovery and exclusion limits set by a 1 TeV LC, even for the most 
challenging channel 777 . While the limits set by a 2 or 3 TeV LC are much higher, and 
the required amounts of data for matching these bounds are too large to be a reasonable 
projection for the LHC reach. 

4 Summary and Conclusions 

In a short summary, we calculate the neutral gauge boson production processes 777, 77Z, 
'-fZZ and ZZZ in the SM and LED model at ILC and CLIC. We investigate the integrated 
cross sections, the distributions of some kinematic variables Myy, and y v . The 5a 
discovery and 3a exclusion ranges for the LED parameters Ms are obtained and compared 
between different channels. It turns out that the effects of the virtual KK graviton enhance 
the total cross sections and differential distributions of kinematical observables generally. 
Among the four processes we considered, e~e + — > ^ZZ or e~e + — > 777 process has the 
most significant LED effect with relatively small or large Ms, respectively. While e~e + — > 
'-ffZ has the least contribution from LED diagrams. With the development of linear 
colliders, more information related to LED effects can be obtained experimentally through 
such important productions. At the 3 TeV CLIC, it is expected that jZZ production 
can be used to explore a range of Ms values up to 7.3-9.3 TeV depending on the number 
of extra dimensions. Through 777 production, we can extend this search up to 9.9 TeV, 
While using 77Z and ZZZ productions, lower bounds on Ms can be found, which are 
6.3 TeV and 7.4 TeV, respectively. 
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